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In a recent Letter [1], Schultz et al. proposed the 
existence of heterofullerenes C n -xN x (40 < n < 50) 
based on the pattern of the peaks in the mass spectrum 
of carbon-nitride nanostructures produced by the arc- 
discharge technique. Using PM3 semiempirical calcula- 
tions, they predicted the most stable structures of the 
N-substituted fullerenes, C33N9 and C44N 6 , to have C3 
symmetry. While some level of symmetry seems to max- 
imize stability, the highest symmetry often sacrifices it, 
as in the case of fullerene Cs4 [2] . Using PM3 and density 
functional theory (DFT), we argue that a highly system- 
atic investigation is needed to determine the most stable 
structures of both pristine and N-substituted fullerenes. 
We demonstrate this explicitly for two particular cases, 
C40 and C44N6, for which the authors of Ref. [1] only 
considered four and two isomers, respectively. 

Mathematically, there are 40 possible isomers for C40 
[3]. We optimized the geometries of the 37 isomers 
with non-vanishing HOMO-LUMO gaps using the PM3 
Hamiltonian and the B3LYP functional with STO-3G, 
3- 21G and 6-31G* basis sets as implemented in Gaus- 
sian 98 [4] . The closed-shell formalism was used and the 
topological symmetry was maintained throughout. Iso- 
mers 1, 30 and 40 (note: the isomer index follows Ref. 
[3]) have zero HOMO-LUMO gaps at the Huckel level of 
theory. Geometric convergence could not be achieved for 
these cases at the PM3 level and no further calculation 
was performed. 

The most stable isomer of C40 is predicted to be iso- 
mer 38 (D2) by all levels of theory. Analysis of the heats 
of formation (not shown) indicated that isomer 38 is not 
among the four isomers studied in Ref. [1], although the 
authors did not identify the isomers. Isomers with low 
relative energies include isomer 39 (Dsd, 10.9 kcal/mol) 
and isomer 31 (C s , 16.1 kcal/mol), see Fig. 1. All other 
isomers have energies higher than 20 kcal/mol. These re- 
sults are in agreement with earlier predictions [5,6]. The 
stabilities from different levels of theory agree with each 
other. 

In Ref. [1], Schultz et al. showed an isomer with C3 
symmetry as the most stable structure for C50 substi- 
tuted with six nitrogen atoms (Fig. 4b in Ref. [1]). Here 
we present a complete search of the Ng-substituted iso- 
mers within the frame of Cso:270 (-D3). Dividing the six 



N atoms into two groups of three atoms and substituting 
different groups of carbon atoms in Cso:270 gives 64 pos- 
sible structures. We first performed calculations at PM3 
and B3LYP/STO-3G levels, which gave similar relative 
energies. Three structures have relative energies lower 
than 16 kcal/mol, whereas all others have relative ener- 
gies higher than 40 kcal/mol. Further calculations at the 
B3LYP/6-31G* level gave relative energies as 0.00 (D a ), 
6.78 (C3) and 9.63 (-D3) kcal/mol, confirming results from 
lower levels. The HOMO-LUMO gaps, 2.70 eV, 2.85 eV 
and 2.58 eV, respectively, are large by DFT standards, 
indicating their kinetic stability. Fig. 2 shows the most 
stable structure of C5o:270 with AVsubstitution having 
D3 symmetry. 

In summary, quantum chemical calculations at both 
DFT and PM3 levels indicate that either a complete sur- 
vey or a more rational selection scheme is needed to de- 
termine the ground state structures of fullerenes. 
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Caption of Figures 

Fig.l: Relative energy of isomers of fullerene C40. The 
isomer index follows Ref. [3] . 

Fig.2: The most stable structure of C 50 :270 with N 6 - 
substitution has D3 symmetry. 
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